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Abstract: We study the twisted Eguchi-Kawai (TEK) reduction procedure for 
large N unitary matrix lattice models. In particular, we consider the case of two- 
dimensional principal chiral models, and use numerical Monte Carlo (MC) simula- 
tions to check the conjectured equivalence of TEK reduced model and standard lattice 
model in the large- N limit. The MC results are compared with the large- N limit of 
lattice principal chiral models to verify the supposed equivalence. The consistency 
of the TEK reduction procedure is verified in the strong-coupling region, i.e. for 
(3 < f3 c where j3 c is the location of the large- N phase transition. On the other hand, 
in the weak-coupling regime (3 > /3 C , relevant for the continuum limit, our MC results 
do not support the equivalence of the large- N limits of the lattice chiral model and 
the corresponding TEK reduction. The implications for the correspondence between 
TEK model and noncommutative field theory are also discussed. 
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1. Introduction 

A renewed interest in reduced matrix models has been motivated by the recent dis- 
covery of their relevance to noncommutative field theories and to superstring matrix 
models, see, e.g., Refs. [[1], 0, f|, ||, ||, |7], ||, |9| and references therein. Reduced matrix 
models were introduced for the study of the large- N limit of SU(N) gauge theories, 
see, e.g., Ref. for a review. Eguchi and Kawai fill pointed out that, as a conse- 



quence of the large- N factorization, one can construct one-point theories equivalent 
to lattice SU(N) gauge theories in the limit N — > oo. The original proposal was just 
the replacement of all link variables U^x) of the Wilson formulation of lattice gauge 
theories with four SU(N) matrices according to the reduction rule U^(x) — > U^. This 
simple reduction procedure was shown to fail in the weak-coupling region, because of 



the spontaneous breaking of a symmetry that is crucial for the equivalence [[12|| . Mod- 
ifications of the original Eguchi-Kawai model have been then proposed |T2|, |T^, |TJ]. 
The most promising one is represented by the so-called twisted Eguchi-Kawai (TEK) 
reduction |13|, [15|, which replaces U^{x) — > D{x)U^D{x)\ where D(x) = Yl^V an d 
are appropriate SU (N) matrices obeying the 't Hooft algebra. 
Since the proposal of the application of the TEK reduction prescription [13|, [0J, 
the numerical checks [|14|, |16|, [L7|, |18| of the validity of the procedure have not been, 
to our opinion, completely exhaustive. We believe that, also in view of the renewed 
interest for matrix models, and particularly for TEK reduced versions, it is important 
to perform an accurate test of the conjectured equivalence with the corresponding 
lattice matrix model. This point may eventually be of relevance for the role of TEK 
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models in noncommutative field theories or in superstring matrix (toy) models, where 
the theory should be analyzed in a limit different from the planar one. 

In this paper we study the application of the twisted Eguchi-Kawai (TEK) re- 
duction procedure to two-dimensional principal chiral models on the lattice |13|, [19], 
20 , 2T| . Two-dimensional chiral models and four-dimensional gauge theories present 



interesting analogies, see, e.g., Refs. |22], |23|, [24], |2^] . For example, they are both 
asymptotically free unitary matrix theories, and their large- N limit is represented 
by a formal sum of planar graphs. In the large N limit, one can apply the TEK 
reduction procedure to both lattice theories. In both cases the validity of the pre- 
scription is demonstrated by showing that the Schwinger- Dyson equations for the 
reduced and the non-reduced theory become formally the same in the limit N = oo. 
In two-dimensional lattice chiral models a stringent test of the TEK reduction pro- 
cedure should be simplified, for essentially three reasons: a better knowledge of the 
large- iV limit of the lattice chiral model, an easier way to define observables in spin 
models, and a simpler reduced model constituted by only one unitary matrix. 

It has been recently realized that TEK reduced model may provide a nonper- 
turbative formulation of noncommutative field theories f|, |5|, || ||, |j. In particular, 
one may establish a correspondence between the TEK reduced model of the lattice 
principal chiral theory and a noncommutative U(l) principal chiral model [26|. This 
is formally achieved in a continuum limit that keeps the noncommutativity param- 
eter 0^ fixed ||, ||, H H) B' wnere @iw scales as a 2 N and a is the lattice spacing. 
Assuming the equivalence of the TEK reduced and the commutative principal chiral 
model in the planar limit, one may rewrite this condition as a double scaling limit 
involving the lattice coupling and iV Therefore, our study of TEK reduced 

model in the planar limit is also useful to clarify the correspondence between TEK 
model and noncommutative field theory, and therefore the possibility of achieving a 
nonperturbative lattice formulation of a noncommutative field theory. 

The plan of the paper is as follows. In Sec. |2] we introduce the original field 
theory and its reduced counterpart, and show on what grounds one expects that it 
coincides with the non-reduced theory in the large N limit. In Sec. |3] we present our 
numerical results obtained by performing Monte Carlo simulations. Finally, Sec. § 
contains comments and conclusions. 



2. TEK reduced chiral models 
2.1 TEK reduced models 

Twisted Eguchi-Kawai reduced models are based on the idea that, when iV — > oo, 
the SU(iV) group becomes so large that it may accommodate the full Poincare group 
as a subgroup, and allow us to find representations of the translation and rotation 
operators among the elements of SU(iV). As a consequence, one may reformulate 
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the full theory in terms of a finite number of matrix field variables defined at a single 
space-time site (or on the d links emerging from the site in the case of a lattice gauge 
theory) and of the above-mentioned representations of the translation group. This 



reformulation is called "twisted Eguchi-Kawai" reduced version of the theory [13, 15 

Consider a lattice field theory with fields 4>(x) in the adjoint representation of 
SU(N). The twisted reduction prescription reads: 

<f>(x) -> D(x) (f) D\x), (2.1) 

where 

D(x) = JIW. (2-2) 

and T M are SU(N) matrices obeying the so called t'Hooft algebra: 

r^r^ — ZpuTvT m (2-3) 

where Z^ v is an element of the center of the group, Z N , 

f2ni \ . . 

Z^ v = exp I jj-n^ I (2.4) 

and n^ u is an integer valued antisymmetric d X d matrix (in d dimensions). is the 
matrix implementing the translation by one lattice spacing in the fi direction. The 
reduced action is simply obtained by substituing (|2.1|) into the action of the original 
lattice field theory: 

Stek^iU^) = -Ls[D(x)<f>D^x)], (2.5) 
and the partition function reads: 

Ztek = J Vcj) exp(-S T EK), (2.6) 

for a fixed value of Z^ u . The expectation value of any functional of the reduced field 
<j) is given by: 



{0(<P))tek = ? / Vcf) 0(0) exp(-S TEK ) (2.7) 

^TEK J 

The correspondence between correlation functions of the reduced model and of the 
original field theory is as follows. Let f(x) be any invariant functional of the field 
4>(x). Then 

field theory = (f[D(x)4>D^(x)]) TE K (2.8) 

The actual check of validity of the reduction procedure is based on the com- 
parison of the Schwinger-Dyson equations of the reduced and original models. This 
procedure however requires some attention, since the limit of infinitely many degrees 
of freedom within the group itself allows the possibility of spontaneous breakdown of 
some of the symmetries which would be preserved for any finite value of N. More- 
over, large iV is a thermodynamical limit: iV must go to infinity before any other 
limit is considered, and sometimes the limiting procedures do not commute. 
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2.2 Principal chiral models 

The two-dimensional principal chiral field theory is defined through the following 



action 



S = ^Jd 2 xTr [d,JJ{x)dJj\x 



[2.9) 



where U(x) are unitary matrix variables. The action is naturally invariant under the 
transformations U — ► VlU and U — > UVr where Vl, Vr G U(N). As amply discussed 



in the literature (see, e.g., Refs. [[23|, |24|, and references therein), two-dimensional 
chiral models and four-dimensional lattice gauge theories manifest deep analogies, 
such asymptotic freedom, large- N planar limit, etc... 

A corresponding lattice principal chiral model can be obtained by the usual 
substitution of the derivative with a finite difference, obtaining 



S 



x fi=l,2 



Tr 



u x ul +lx + u x+u ut 



(2.10) 



where (3N = 1/T. In view of a large- iV analysis one may consider both SU(N) and 
U(N) models, since they are expected to reproduce the same statistical theory in 
the limit iV — > oo, which is meant at fixed (3. Moreover, SU(N) and U(N) chiral 
models have the same continuum limit at any finite N > 2. 



For a recent review of results concerning chiral models see, e.g., Ref. ||24|| . Here we 
only mention the presence of a peculiar large- N phase transition |^7|, [3(J which 

manifests itself with a power-law singularity in the N —> oo limit of the specific heat, 
i.e. Coo ~ \/3 c ~ (3\~ a , and with a peak in the specific-heat at finite iV that becomes 
sharper and sharper with increasing iV |28|, [SIJ. The analysis of the N = oo strong- 
coupling series [29j provides the estimates f3 c = 0.3060(4) and a = 0.27(3), while the 
extrapolations to iV = oo of the position of the specific-heat peak, determined by 
Monte Carlo simulations for SU(N) and U(N), gives (3 C = 0.3957(3). 

2.3 The reduced version: symmetries and equivalence 

Implementing fl2.1p into Q2.10D gives the TEK reduced lattice chiral model action: 



Stek = -PN Tr [ Ur » U K + h - c - 



H=l,2 



In the present case, of equation ( |2.4| ) is generically, in d — 2, of the form 

M' 



—M 



M e z. 



(2.11) 



(2.12) 



For a given iV and M the solution to (f3T[) is provided, up to global SU(N) trans- 
formations, by the N x N shift and clock matrices 



q(M) 



5 



C 



i+M,j i 



e ¥H) 5 .. 



(2.13) 
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The two matrices T M are given respectively by S and C. 

While the original field theory is invariant under the SU(N)l x SU(N)r group, 
the reduced theory is invariant under the following two symmetries: 



U -> z ■ U z e Z. 
U -> D{x) U D\x). 



(2.14) 
(2.15) 



The first symmetry corresponds to the contraction of the original symmetry of the 
model to the center of the algebra of the symmetry group, while the second one is 
reminiscent of the translational invariance of the original theory. 

In order to show the equivalence of the reduced theory with the original field 
theory in the large N limit, the starting point is given by the Schwinger-Dyson 
equations, i.e. the equation obtained from a generic correlation function in the 
original theory (the reduced case is the obvious translation using Eq. Q2.1| )) 



by the following steps: 

1. substitute with X a U^, where A a is a SU(N) generator; 

2. implement the change of variables U — *■ (1 + ia\ b )U in the integral (implicitly 
taken in the average), using the Haar measure invariance; 

3. contract the a, b indices and extract the linear part in the parameter a. 

By comparing the two sets of equations obtained in the original and in the reduced 
theory, one finds manifestly corresponding terms and also some extra terms in the 
reduced theory. These extra terms exactly cancel if the symmetry U — ► D{a)UD\a) 
of the reduced model is not spontaneously broken fl9|| . Therefore, the two models 
possess the same Schwinger-Dyson equations if the above mentioned symmetry is 
not broken in any regime, and in particular in the continuum limit. The point is 
then whether this implies the equivalence of the two models, i.e. that they have 
the same expectation values of invariant quantities. Since the large N equivalence 
of the twisted reduced models is based on the above arguments even for the lattice 
gauge theory case, it is an important task to verify, if the invoked symmetry holds, 
whether the Schwinger-Dyson equation based demonstration of the equivalence is 
indeed conclusive and sufficient or not. 

The equivalence of the lattice chiral model and the reduced TEK model in the 
N = oo limit is not actually questioned in the strong-coupling region, i.e. for suffi- 
ciently small values of (3. In this region the symmetries of the TEK reduced model are 
unbroken and therefore the equivalence of Schwinger-Dyson equations holds. One 




(2.16) 
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may then argue that this is sufficient to yield the equivalence of the reduced and 
nonreduced models. In the strong-coupling domain, i.e. within the convergence ra- 
dius of the strong-coupling expansion, expectation values are analytic in the coupling 
P, and their boundary value at (3 = can be easily computed. As a consequence, 
it is formally possible to solve the Schwinger-Dyson equations in terms of strong- 
coupling series by iterative methods. Since the N = oo boundary values at (3 = 
of the lattice chiral model and the reduced TEK model coincide, equal Schwinger- 
Dyson equations imply the formal equivalence within the convergence radius of the 
strong coupling expansion 1 . Problems may occur in the weak-coupling side of a 
large- N phase transition. 



It has been also argued that |H| the reduced TEK model is, in a sense, equivalent 
to a theory on a box of size L = N. In the large- N limit the finite- N corrections 
should beO(l/iV 2 ), just as in the SU(N) lattice chiral theory. Since N 2 = L 2 , finite- 
iV corrections can be seen as finite volume corrections. Therefore, one may expect 
that the asymptotic large- N regime is reached for N >> £, where £ is the N = oo 
correlation length at the given value of (3. Note that in TEK reduced models the 
large-iV and thermodynamic limits are connected and approached simultaneously, 
while in the standard lattice models one first performs the L — > oo limit and then 
the large- N limit. In a general field theory there are no a priori reasons why these 
two limits should commute. 



3. Numerical results 

In order to study the large-iV limit of the TEK model, we have performed Monte 
Carlo simulations for several values of N, up to iV = 200, and for values of (3 in 
the strong and weak coupling regime. We used a Metropolis algorithm to update 
the SU(N) matrix U. Trial matrices were selected by multiplying the actual matrix 
U by a random SU(2) matrix embedded in SU(N), choosing randomly among the 
N(N — 1)/2 577(2) subgroups. More precisely, once the £77(2) subgroup is randomly 
chosen, we performed ten Metropolis hits with an approximate acceptance of 50%. 
Each SU(2) updating requires 0(N) operations. The number of £[/(2)-subgroup 
updatings per run was O(10 9 ). We should also mention that for large values of N, 
N > 60 say, we observed some problem of thermalization when using completely ran- 
dom configurations (for example, constructed by multiplying N(N — 1)/2 completely 
random SU(2) matrix embedded in a N x N and associated to different subgroups) 
as starting point of our simulations. This occurred in both the strong- and weak- 
coupling regions, worsening with increasing (3. So, as starting point we used either 
moderately random matrices or the unity matrix. It is well known that a simple 

1 Notice that the convergence radius of the strong coupling expansion does not necessarily coincide 
with the large-TV transition point /3 C , but it may be smaller. 
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Metropolis algorithm does not provide a particularly efficient method to simulate 
a statistical system. Our choice was essentially due to the fact that the reduced 
action is quadratic in the matrix variable U. Moreover, it does not lend itself to a 
linearization by introducing new matrix variables, as in the case of the reduced TEK 
gauge theory fTTjj . 

In our simulations we measured the internal 
energy 

E = l-±(TrUr^rl), (3.1) 

the quantity corresponding to the two-point func- 
tion G(x) = (Tr U(0)U(xY) using the equiva- 
lence (|2.8| ), and in particular its low- momentum 
components, such as the (zero-momentum) mag- 
netic susceptibility 

X =\TrU\ 2 . (3.2) 

In order to check for possible spontaneous break- 



AT 


TP 


X 


10 


0.6155(2) 


12.30(2) 


16 


0.6346(1) 


12.47(2) 


20 


0.6409(1) 


11.46(4) 


30 


0.64766(7) 


9.49(3) 


40 


0.65046(8) 


8.44(3) 


50 


0.65159(5) 


7.94(3) 


60 


0.65226(4) 


7.65(3) 


80 


0.65285(3) 


7.34(4) 


100 


0.65329(4) 


7.21(4) 


150 


0.65349(4) 


7.08(5) 



Table 1: Internal energy E and 
magnetic susceptibility \ f° r P = 
0.28. 



ing of the symmetries (|2.14j) and ( |2.15|) , we con- 
sidered the two quantities 



Si 



N 



(TrlO, 



S 2 =l(Tr£TWt r t). 



(3.3) 
(3.4) 



Indeed, Si is not invariant under the transformation fl2.14|) (but it is invariant with 
respect to ( |2.15|) ), while S2 is not invariant under the transformation (|2.15| ) (but it 
is invariant with respect to ( |2.14|) ). 

Information on the large- N limit of lattice chiral models in the strong coupling 
region can be obtained from the N = 00 strong-coupling expansion and extrapolation 
of Monte Carlo results for the SU(N) and U(N) unitary groups. N = 00 strong- 
coupling series are reported in Refs. [|22], |27|, |28|, [29|]. For example, the internal energy 
is known to 17th order and the magnetic susceptibility to 15th order. Monte Carlo 
simulations for relatively large values of N (up to N = 30) can be found in Refs. 
[^, 3(J. It is useful to consider both SU (N) and U (N) data because their expectation 
values aproach the N = 00 values from opposite sides, and allow us to have a better 
control of the large- N limit. 

As already mentioned, the equivalence of reduced and nonreduced models is not 
questioned in the strong-coupling region. However, Monte Carlo simulations in the 
strong coupling region are still useful to study the convergence of the large- N limit 
of TEK, which is expected to be controlled by powers of N~ 2 . 
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After verifying a general agreement between the large- N TEK results and the 
N = oo strong-coupling series for small (3, we have performed longer runs for (3 = 0.28 
in order to obtain an accurate test. This value of (3 is not too far from the large- N 
transition point (3 C ~ 0.306. The large- N limit of E and x of the lattice chiral model 



can be estimated from the analysis of the corresponding strong-coupling series [2E 
obtaining 2 = 0.65370(4) and Xoo = 7.02(1). Extrapolations of the Monte Carlo 
results reported in Refs. ^U|, using the fact that the corrections are expected 
to be powers of N~ 2 , provide consistent results, i.e. 3 Xoo = 7.00(2) for SU(N), 
Xoo = 7.02(2) and E^ = 0.65373(4) for U(N). Moreover, we mention that the large- 
N limit of the second-moment correlation length is ^ ~ 1.53. In Table |l] we report 
E and x as obtained by our simulations for several values of N, up to iV = 150. In 
Figs. [L]and|2|we plot them versus l/N 2 , which is the order of the expected asymptotic 
corrections. There we also show the MC data obtained from simulations of SU(N) 
and U(N) lattice chiral models. The figures show clearly that the correct N = oo 
limit is approached. To be more quantitative, fitting the data for iV > 30 allowing 
for a 0(N~ 2 ) correction, we obtained 

E « 0.65377(3) - (3.5) 

x „ 7.02(2) + ^f^, ,3.0) 

with acceptable x 2 /d-o.f ~ 1. This analysis shows that the correct large- N limit 
is approached with corrections of 0(N~ 2 ), as expected. Note, however, that the 
0(1/N 2 ) corrections are much larger than those found for the lattice chiral models, 
by about a factor 10 2 for x, see the footnote 

As already mentioned, lattice unitary chiral models present a peculiar large- N 
phase transition, which is characterized by a power-law singularity of the specific 
heat at iV = oo |2£|. At finite N, numerical simulations show that the peak in the 
specific-heat becomes sharper and sharper with increasing iV |3(J. A peak in the 



specific heat is also observed in the TEK reduced model: it becomes sharper and 
sharper with increasing N. The iV — > oo extrapolation of its location provides the 
estimate (3 C ~ 0.306, in agreement with the results for the lattice chiral models. 

These results provide an accurate check that for (3 < j3 c the TEK model shows 
the expected behavior, i.e. it converges to the same large-iV limit of lattice chiral 
model, with corrections of 0(N~ 2 ). We should also observe that, due to the large 

2 Following Ref. p9| , we used Dlog-Pade approximants to resum the series of dE/dfl ~ (/3 C — ft)~ a 
and dlnx/dp ~ (j3 c — (3)~ a with j3 c « 0.306 and a w 0.27. The internal energy and the magnetic 
suceptibility are then obtained by integrating with respect to j3. 

3 For example, the simplest extrapolation using a straight line a + bN~ 2 passing by the data 
for the largest values of N (i.e. x(N = 30) = 7.032(5), x(N = 21) = 7.069(5) for SU(N) and 
X (N = 21) = 6.972(8), X (N = 15) = 6.924(7) for U(N)) gives xW = 7.00(1) + 32(6)iV- 2 for 
SU(N) and x(N) = 7.02(2) - 22(5)A^" 2 for U(N). 
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Figure 1: The internal energy versus N~ 2 for (3 = 0.28. We also show the large- iV limit 
.Boo as obtained from the analysis of the N = oo strong-coupling series, and results from 
simulations of SU(N) and U(N) lattice chiral models. 



corrections, the convergence to the N = oo limit is rather slow, much slower than 
the convergence of the lattice chiral models. Comparing with the numerical effort 
to reach a similar accuracy in simulations of lattice chiral models, the reduced TEK 
model does not appear to be convenient to numerically investigate the large- N limit 
of chiral models, at least in the strong-coupling domain. 

Let us now turn to the weak-coupling region, i.e. for (3 > (3 C . In this case 
information on the large- N limit of lattice chiral models comes only from the MC 
simulations of Refs. 1 28, |30 ] . We performed runs for (3 = 0.31 and (3 = 0.32, for which 
the N = oo correlation length are estimated to be £ ~ 4 and £ ~ 6, respectively. MC 
simulations up to iV = 200 for these values of f3 should provide a convincing check 
of the large- N convergence, which should be controlled by the ratio iV/£ according 
to the arguments mentioned in Sec. |2.3| . The results of the simulations for E, Xi 
and the two quantities S\ and S2, cf. Eqs. ( |3T3"| ) and (|3.4j), are reported in Table 0. 
They should be compared with the the extrapolations to iV — > 00 of the MC results 
reported in Ref. [|^, [ftj, that are Xoo = 34.1(2) and E^ w 0.519 for (3 = 0.31, and 
Xoo ~ 65 and E^ w 0.485 for (3 = 0.32. Our MC data are also plotted in Figs. | and 
|], where we also show the MC results for SU(N) and U(N) lattice chiral models, 



taken from Refs. |28| , pO]] . While the data for the energy seem somehow to approach 
the corresponding E^ (although it is not evident that the corrections are 0(N~ 2 )), 
those for \ do not show any evidence of convergence to Xoo- X appears to increase 
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Figure 2: The magnetic susceptibility versus 

AT- 2 for (3 = 0.28. We also show the large- TV 
limit Xoo as obtained from the analysis of the N = oo strong-coupling series, and results 
from simulations of SU(N) and U(N) lattice chiral models. 



with N up to iV = 200. Similar results are obtained for (3 = 0.32. We have also 
verified that the reduced two-point function does not exponentially decrease at large 
distance, explaining the behavior of its space integral x- Finally, note that the results 
for Si and S 2 , see Table ^ are always consistent with zero, so that they do not give 



any indication of breaking of the symmetries (|2.14|) and (|2.15| ). 

Our results apparently contradict earlier numerical works [J21], |2"H] , which were 
based on a much smaller statistic, and where the comparison was essentially restricted 
to the internal energy and the quantity Si. 



4. Conclusions 

We have investigated the application of the twisted Eguchi-Kawai reduction proce- 
dure to two-dimensional principal chiral models on the lattice, as general testcase of 
the TEK reduction procedure. 

In the strong coupling region, i.e. for (3 < /3 C , the large-iV limit of the reduced 
theory appears to reproduce the results available for the original field theory in 
the large- N limit, and that the TEK reduced model approaches the large- N limit 
with power corrections of iV~ 2 . On the other hand, the validity of the reduction 
prescription apparently breaks down for (3 > (3 C . In particular, it is not supported 
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N 


E 


X 


Si 


Re S 2 


0.31 


10 


0.5251(1) 


22.06(2) 


-0.0011(8) + i0.0004(9) 


0.00002(3) 




16 


0.5329(2) 


35.46(6) 


-0.001(2) +i0.001(2) 


0.00000(3) 




20 


0.5342(2) 


43.9(2) 


-0.004(3) + i0.005(3) 


-0.00001(5) 




30 


0.5337(2) 


64.0(2) 


0.003(2) -i0.005(2) 


-0.00002(3) 




40 


0.5317(3) 


83.3(6) 


0.001(4) + z0.005(3) 


0.00006(4) 




60 


0.5265(3) 


123(2) 


-0.003(8) + i0.003(6) 


0.00011(5) 




80 


0.5236(2) 


157(2) 


0.007(6) -i0.001(5) 


0.00003(4) 




100 


0.5216(2) 


188(3) 


-0.013(9) -i0.014(8) 


0.00005(5) 




150 


0.5204(1) 


211(10) 


-0.011(7) +i0.004(10) 


0.00002(6) 




200 


0.5200(1) 


223(4) 


-0.010(9) + i0.004(9) 


0.00004(4) 


0.32 


40 


0.4902(1) 


148.2(7) 


-0.009(6) -i0.002(5) 


0.00001(4) 




60 


0.4880(2) 


235(4) 


0.03(2) + i0.03(2) 


-0.00003(7) 




80 


0.4869(2) 


331(3) 


-0.01(2) - i0.03(2) 


-0.00011(7) 




100 


0.48656(4) 


404(2) 


0.03(2) - i0.02(2) 


0.00001(4) 




150 


0.48629(8) 


565(7) 


0.01(2) - i0.06(2) 


0.00013(8) 



Table 2: We report the intenal energy E, the magnetic susceptibility x> and the quantities 
S\ and 52, defined in Eqs. (|3^ ) and (|Q|), for = 0.31,0.32 and several values of N. 



by the MC results for the magnetic susceptibility, that do not show evidence of 
convergence up to iV = 200 and for (3 values corresponding to large- iV correlation 
lengths £ ~ 4 and £ ~ 6 (estimated extrapolating MC data for lattice chiral models) . 
Of course, we cannot exclude that the correct asymptotic regime is eventually set for 
larger values of N, N ^> 10 2 say. Although we note that already for smaller values of 
(3 < j3 c the approach to the large- limit is rather slow, in order to explain our results 
we should suppose a strong slowing down of the convergence with increasing (3. In 
this regard, we should also add that the equivalence of the theories in the continuum 
limit requires a uniform convergence to A^ = oo when approaching the continuum 
limit. Another possibility is that our numerical simulations have been trapped by 
nontrivial extrema of the action, which has nothing to do with the correct limit 
related to the chiral field theory. However, our numerical checks definitely exclude 
it at least for A^ < 60. 

We have also checked for possible spontaneous breakings of the symmetries ( |2.14j ) 
and Q2.15 ) of the TEK reduced model, by computing the quantities 5*1 and 5*2, cf. 
Eqs. (|3.3j ) and (|3.4| ), which are respectively not invariant with respect to the sym- 
metry (|2.14 ) and (|2.15|) . Their expectation values appear to vanish, thus suggesting 
that the symmetries of TEK are not spontaneously broken, although one cannot 
exclude breaking of symmetries associated with more complicated order parameters. 
We recall that one of such symmetries, i.e. the symmetry ( p. 15] ), is essential for 
the demostration of the equivalence of the N = oo Schwinger- Dyson equations. As- 
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Figure 3: The internal energy versus N~ 2 for = 0.31. We also show the MC data for 
SU(N) lattice chiral models, and their large-iV extrapolation. 



suming that the Schwinger- Dyson equations become formally equal for N = oo, the 
question to be addressed is whether this is sufficient to identify a theory. Our numer- 
ical results suggest that this could not be the case, under the hypothesis that one can 
recover the N = oo theory from the finite N behaviour of the reduced theory in its 
N —* oo limit. In other words, the N = oo TEK reduced theory may be still equiva- 
lent to the large- N limit of lattice models, but the limit is singular, so that it cannot 
be reached from finite values of N. This issue is of relevance in any application of 
the TEK reduction procedure, as in the case of the lattice gauge theory. 

The apparent failure of the TEK reduced theory to reproduce the planar limit of 
lattice chiral models does not necessarily prejudice the possibility of considering the 
TEK theory system in a different limit, possibly relevant to the fields where twisted 
Eguchi-Kawai reduced models were recently applied: noncommutative field theories 
and superstring matrix models. Indeed, TEK reduced chiral models with symmetry 
group U(N) are formally equivalent to a U(l) noncommutative theory f26|. This 
equivalence requires a nontrivial continuum limit that is different from the planar 
limit we have considered in this work, since one needs to keep the noncommutativity 
parameter 6^ fixed, and 6^ scales as a 2 N where a is the lattice spacing. Thus, the 
relevant continuum limit should be realized for a — > and N — > oo keeping a 2 N 
fixed. The lattice spacing is somehow controlled by the coupling (3, which should 
be appropriately tuned to realize the limit a — > 0. The relation between a and (3 is 
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Figure 4: The reduced magnetic susceptibility for (3 = 0.31 versus N~ 2 . We also show 
the MC data obtained by simulating SU(N) and U(N) lattice chiral models; the lines 
represent their large-iV extrapolations. 



therefore essential to define the continuum limit associated with the noncommutative 
field theory. However, the apparent failure of TEK model to reproduce the planar 
limit of the principal chiral theory does not allow us to use the known asymptotic 
behavior of the commutative theory in the continuum limit, for which asymptotic 
freedom and renormalization group imply a ~ exp(— 871 P) in the continuum limit 
a — > and for N = 00. Of course, this point deserves further investigation. 
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